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Abstract 

Bound state properties of the negatively charged Ps _ ion (or e _ e + e _ ) are discussed. The 
expectation values of operators which correspond to these properties have been determined with 
the use of the highly accurate wave functions constructed for this ion. Our best variational energy 
obtained for the Ps _ ion is E = -0.2620050 7023298 0107770 40051 a.u. Annihilation of the 
electron-positron pair(s) in the negatively charged Ps _ ion (or e _ e + e _ ) is considered in detail. By 
using accurate values for a number annihilation rates r n7 , where n = 1, 2, 3, 4 and 5, we evaluated 
the half-life r a of the Ps _ ion against positron annihilation (r a = ~ 4.793584140 • 10 -10 sec). 

Photodetachment of the Ps _ ion is considered in the long-range, asymptotic approximation. The 
overall accuracy of our photodetachment cross-section is very good for such a simple approximation. 
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I. INTRODUCTION 


The main goal of this short communication is to perform computational and theoretical 
analysis of the annihilation of electron-positron pair(s) in the negatively charged Ps _ ion. 
Another aim is to evaluate the photodetachment cross-section of the Ps _ ion in the long- 
range, asymptotic approximation. Our analysis is based on some recent results of highly 
accurate computations performed for the ground (bound) l 1 ^—state in the Ps _ ion, which 
is also designated as the e _ e + e _ ion, or e + e^ ion. Stability of this three-body system 


was predicted by Ruark [1], First variational calculations of the ground state in the Ps“ 
ion were performed by Hylleraas in 1947 j^. This ion is of interest in various branches of 

physics of high-temperature 


physics, including solid state physics Q , astrophysics fl.fl.B, 
plasmas, etc. Note that the Ps _ ion has been created experimentally by Mills almost forty 
years ago 7]. Review of the most recent experiments performed for the Ps~ ion can be 
found in jj and |9] which also contain a large number of useful references. 

As is well known from the general theory of bound states in the Coulomb three-body sys¬ 
tems with unit charges jioj] this ion has only one stable state, which is the ground HA—state, 
or 1 1 S'(L = 0)-state. In general, to describe the bound IPS'—state in this ion we can restrict 
ourselves to the non-relativistic LS— approximation, since all lowest-order relativistic and 
QED corrections are very small for this ion. The non-relativistic approximation means that 
the wave function T can be determined as the solution of the non-relativistic Schrodinger 
equation H4> = AT, where E < 0, for three-particle quasi-atomic (or Coulomb) systems. 
The non-relativistic Hamiltonian H of the Ps” ion takes the form 


H = 


2 m P 




e 

t-m 


e 

r- S2 


e 

T 21 


( 1 ) 


where h — is the reduced Planck constant (also called Dirac constant) and m e is the 
electron mass and —e is the electric charge of an electron. In this equation and everywhere 
below in this study the subscripts 1 and 2 designate two electrons (e _ ), while the subscript 
3 always denotes the positron (e + ) with the mass m e (the same electron mass) and positive 
electric charge +e, or e. In addition to the ‘numerical’ indexes (1, 2, 3) in some cases 
we shall designate electrons by using the notaiton while the notation ‘+’ always means 
positively charged positron. In Eq. (JTj) the notations =| r* — |= -r^ stand for three 

interparticle distances (= relative coordinates) which are the absolute values of differences 
of the Cartesian coordinates of the three particles. Note that each relative coordinate ry,- 
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is a scalar which is rotationally and translationally invariant. However, these coordinates 
are not truly independent, since e.g., | r 32 — r 3 i |< r 2 1 < r 32 + r 3 i. This produces a number 
of problems for computations of three-particle integrals in these coordinates. To simplify 
such calculations it is better to apply a set of three perimetric coordinates u 3 ,u 2 ,u 3 which 
are simply related to the relative coordinates: tq = | (r^ + r^. — r t] ), while inverse relations 
take the form r y - = iq + Uj. Three perimetric coordinates u\,u 2 ,u 3 are independent of each 
other and each of them varies between 0 and + 00 . The Jacobian of the transition u 7 

D UltU2 , U3 (r 32 , r 3 i, r 2 i) is a constant which equals 2. 

Note also that in this study only atomic units h = 1, | e |= 1 ,m e = 1 are employed. In 
these units the explicit form of the Hamiltonian iJ, Eq. (JTJ), is simplified to the form 

1 r 


H = — ^ |V^ + V2 + V3 


111 

r 3 1 r 32 r 21 


( 2 ) 


Note that the Hamiltonian, Eq.([2]), does not contain any ratio of masses and/or electric 
charges. It follows from here that the Ps _ ion plays a central role in the general theory of 
Coulomb three-body systems with unit charges (for more details, see, e.g., 10]). 

To solve the non-relativistic Schrodinger equation ET = ET for the Ps _ ion, where 
E < 0, and obtain highly accurate wave function(s) we approximate the unknown exact 
solution of the non-relativistic Schrodinger equation with some efficient and fast convergent 
variational expansions. The best of such expansions is the exponential variational expansion 
in the relative coordinates r 32 ,r 3 i,r 2 i, or perimetric coordinates Ui,u 2 ,u 3 . For the ground 
state of the Ps _ ion the explicit form of this expansion is 


1 


N 


V I ; = -(1 + Pi 2 ) C i ex p(—cqr 32 - Apn - 7* r 2i) 

2 i=i 

1 , N 

= -(1 + P12) Ci ex p[-(«i + 3 - {oti + 7 i)u 2 - (A + 7 i)w 3 ) 

2 i=i 


( 3 ) 


where the notation P\ 2 stands for the permutation operator of identical particles, Ci (i = 
1,2,..., N) are the linear parameters of the exponential expansion, Eq.Q, while cq, A and 
7 i are the non-linear parameters of this expansion. The non-linear parameters must be varied 
in calculations to increase the overall efficiency and accuracy of the method. The best-to- 
date optimization strategy for these non-linear parameters was described in llj, while its 
modified version is presented in fl2 |. The 3iV— conditions cq + A > 0 ,cq + 7 * > 0 , A + 7 « > 0 
for i = 1,2, ...,1V must be obeyed to guarantee covergence of all three-particle integrlas 
needed in computations. 
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II. EXPECTATION VALUES 


By using the highly accurate, variational wave function T constructed for the ground 
ids'—state of the Ps _ ion we can determine the expectation value of an arbitrary, in principle, 
self-adjoint operator X. This is written in the following general form 

(T | X | T) 


(*) = 


(4) 


(T | T) 

Formally, without loss of generality below we shall assume that our wave function has a unit 
norm, i.e. (T | T) = 1 (see, discussion in jl31 ) . The total energy is the expectation value 
of the Hamiltonian H , Eq. (JTJ) , i.e. E = (T | H | T). The total energies E of the ground 
l 1 S—state of the Ps _ ion determined for different trial wave functions can be found in Table I. 
Other possible choices of operators X in Eq.() lead to the different bound state properties, 
or properties, for short. A number of bound state properties were determined in earlier 
computations of the Ps~ ion (see, e.g., [bfj], (eJ, jit]] and references therein). In this study 
we present a large number of bound state properties of the Ps _ ion determined to very high 
numerical accuracy (‘essentially exact’). They can be found in Table II (in atomic units). 
Physical meaning of many of these properties is clear from the notations used in Tables I 
and II. For instance, the notation (ry,) stands for the expectation value of the linear distance 
between particles i and j. Another notation (5^) = (<5(rj — ry)) denotes the expectation value 
of the (Dirac) delta-function between particles i and j, while (£ 321 ) = (h(r , — iy)h(iy — r fc )) 
is the expectation value of the triple delta-function. In general, the expectation value of 
each delta-function is the probability to locate two (or three) particles inside of one small 
sphere with the radius R ~ A e = aao, where A e is the Compton wavelength of electron, 


a = jP ~ is the fine structure constant and a 0 is the Bohr radius (see below). 


For the Ps _ ion the expectation value of the electron-positron delta-function _) deter¬ 

mines a number of annihilation rates, including the two- and three-photon annihilation rates 

(see below). The expectation value of the triple delta-function (£ 321 ) = _) is important 

to predict the one-photon annihilation rate Ti 7 . It is clear that some reliable criteria are 
needed to check the overall quality of the computed expectation values of delta-functions. 
In reality, we can introduce such criteria by considering the coincidence of the computed 
and predicted cusp values between each pair of particles^ It was shown in early papers on 
Coulomb systems, including atoms and molecules jl7], [l 8 |, that the following expectation 
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value: 


= (%) = 


<* I A I *> 


(5) 


I ^(r*i) I 

is always finite and its numerical value equals to cuq-j m 'f b , where g l; g 7 - are the corresponding 
electrical charges of particles, while mi,rrij are their masses. The expectation value ry, is 
called the cusp between two Coulomb particles i and j. The coincidence of the computed 
expectation value of the cusp u t j = (%) with its expected value, i.e. with q % q 3 , 
indicates the overall quality of the expectation value of the inter-particle delta-function. 
In actual applications to Coulomb systems this criterion works very well. Our computed 
and expected cusp values are presented in Table II. As one can see from Table II numerical 
coincidence between the predicted and computed expectation values of the electron-positron 
and electron-electron cusps can be considered as very good. The predicted value of the 
electron-positron cusp for the Ps - ion equals -0.5 a.u., while for the electron-electron cusp 
one finds 0.5 a.u. Unfortunately, there is no similar criterion for the triple delta-function $ 321 , 
since the corresponding ‘three-particle cusp’ is infinte for an arbitrary Coulomb system (see, 
e-g., EL El and references therein). However, for trial functions with the finite number 
of regular (or non-singular) basis functions the three-particle cusp z / 123 can be defined (this 
value is finite) and used in computations EL 

The notations T tJ stand for the expectation values of the interparticle cosine-functions 
which are defined traditionally: 


Tij = {cos(r, fe A r,fc)) = {— • —) 


( 6 ) 


Tik Tjk 

where ( i,j,k ) = (1,2,3), while the notation r ifc A r j k denotes the angle between the two 

vectors and r^. For the Ps _ ion there are two independent r t j values, i.e. r 2 1 (or r_) 

and T 31 (= 7 ~ 3 2 , or iq_). These two values are shown in Table II. For an arbitrary three-body 

system the sum of the r 2 i,r 3 i and r 32 values is represented in the form 


7-21 + 7-31 + r 32 = 1 + 4 (/) 


(7) 


where the notation (/) is the following expectation value 

/*+00 /*+oo r 

(/) = / / / 

Jo Jo J 1 


+00 /-+oo /■r 3 2+r 3 i U1U2U3 

^(r 32l r 3 i, r 2 i)-'h(r 32 , r 3 i, r 2 i)r 32 r 3 ir 2 idr 32 dr 3 idr 2 i (8) 

0 ^0 J |f 3 2—f 3 l | ^32 r 31 r 21 

/*+00 /*+00 /*+oo 

= 2 / / / 'fy(u 1 ,U2,U3)uiU2U3^(ui,U2,U 3 )du 1 dU2du 3 

Jo Jo Jo 
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where r 32 ,r 3 i,r 2 i are the three relative coordinates and ui, 112,113 are the three perimetric 
coordinates, respectively. For the symmetric two-electron Ps _ ion the equality, Eg. (IT)) . takes 
the form r 2 1 + 2r 3i = 1 + 4(/). In applications to real systems the equality, Eq.(J7]), can be 
considered as an additional test of the trial wave function. Indeed, the expectation values 
in the both sides of Eq.([7J) can be computed separately. The coincidence of the results in 
both sides of Eq. (]7]) indicates correctness of the procedure used. 

In general, there are a number of additional equalities for the expectation values deter¬ 
mined for an arbitrary three-body system (or few-body system). For instance, consider the 
following identity r 32 = r 31 + r 12 . From this identity one finds r| 2 = + r 21 — 2r 31 • r 21 . 

This leads to the following identity for the expectation value (r 31 • r 21 ) 

( r 3 i • r 2 i) = ^((r 31 ) + (r 21 ) - (r 32 )) (9) 

Two other analogous equalities for the (r 32 • r 12 ) and (r 32 ■ r 31 ) values can be derived by 
using the same approach. Another important example follows from the conservation of the 
momentum in an arbitrary three-body system. In the general case, for an arbitrary wave 
function T written in the relative coordinates we can write (Pi + P 2 + P 3 )^ = 0. It is 
assumed here that the wave function T depends upon three relative coordinates which are 
translationally (and rotationally) invariant. From this equation one Ends three different 
identies between corresponding expectation values, e.g., 

(Pi ' P2> = \((pI) ~ ( pi) - (pi)) (10) 

and 

(Pi • Ps) = \((pI) ~ (p\) ~ (pD) (11) 

For the Ps _ ion the first and second expectation values in the right-hand side of Eq. (fill) are 
equal to each other (particles 1 and 2 are identical electrons). Therefore, in this case one 
finds (pi • p 3 ) = (p 2 • p 3 ) = — (pi) < 0. In atomic units for the Ps _ ion we have p& = —?V/, : 
(k = 1, 2, 3), where 1 is the imaginary unit and V& is the gradient operator for k —th particle 
(also called the Hamilton operator). Additional relations between other expectation values 
can be investigated analogously. 

Note that the current accuracy achieved in calculations of a large number of bound 
state properties of the Ps - ion is very high. This follows from Table II where only stable 
decimal digits are shown for each property. Briefly, we can say that many bound state 
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properties of the Ps~ ion are now known to the level which cannot easily be improved in 
future computations, or such an improvement has no direct physical sense. By using the 
expectation values from Table 11 we can evalauate some fundamental properties of the Ps - 
ion. The main interest for the Ps - ion is related to the positron annihilation, or annihilation 
of the (e~, e + )— pair(s). In reality, we have a number of different few-photon annihilation 
rates, including two- and three-photon annihilation rates. These values are evaluated below. 
Note also that numerical evaluations of the relativistic and lowest-order QED corrections for 
the Ps” ion can be found in y (see also [21,]). It should be mentioned, however, that such 
evaluations were not completed. Moreover, some of the values determined in [l 6 ( must be 
re-calculated to better accuracy and with the use of different algorithms. It is also clear that 
by using our current methods we cannot finish our numerial evalution of the lowest-order 
QED correction for the Ps - ion, since this ion has no infinitely heavy, central nucleus. This 
problem is briefly discussed in the Conclusion. 

III. POSITRON ANNIHILATION 


Now, let us discuss the process of positron annihilation, or annihilation of the electron- 
positron pair(s), in the three-body Ps - ion. As mentioned above this process is of great 
interest in numerous applications, including astrophysics, solid state physics, etc. It is also 
important for workability of many technical devices, e.g., modern communication devices. 
In the general case, annihilation of the electron-positron pair(s) in the Ps - ion proceeds with 
the emission of a number of photons (see, e.g., 22 ], 23 ]), e.g., 


e + + e = 71 + 72 + • • • + 7 k 


( 12 ) 


where 7 *. (k = 1 ,..., K) are the emitted photons and K is the maximal number of such 
photons. Each of the annihilation processes has its unique annihilation width, or annihilation 
rate Tfc 7 . For the Ps - ion the following annihilation rates are important in applications: 
^ 7 ^ 37 ^ 47^57 and Ti 7 (here they are ordered by their numerical values). Formulas for 


these annihilation rates were discussed in a number of earlier studies (see 


, e.g., jpj 


and 


21 ] 


and references therein). For instance, the known analytical expression for the r 27 rate is 
written in the form 


r 27 (Ps ) = rnic^ca, 


4„„ —1 
0 


r a 

/ TT" \ 

1 — 

5-) 

L 7r 

( 4 J\ 


100.34560545419 • 10 9 (<5+_) sec' 1 (13) 
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where a = = 7.2973525698 • 10“ 3 (~ -yAj is the dimensionless fine structure constant, 

c = 2.99792458 ■ 10 8 m ■ sec~ l is the speed of light in vacuum, and the Bohr radius a$ equals 
0.52917721092 • 10 -10 m 25]. Also, in this formula h = A i s the reduced Planck constant (or 
Dirac constant) and n is the total number of electron-positron pairs in the polyelectron e+ e~, 
i.e. n = pq. For the three-body Ps - ion we have n = 2. In Eq. (fT3ji and in formulas below 

the notation (<5_|_) stands for the expectation value of the electron-positron delta-function 

determined (in atomic units) for the l 1 S'-state of the Ps - ion. The formula, Eq. (fTTji . also 


includes the lowest-order radiative correction 


24] (or QED correction in modern language). 


By using this formula for the r 27 rate we can evaluate the four-photon annihilation rate r 47 , 
since the relation between these two values takes the form fsoj ] 


/ft\ 2 „ 

r.i 7 (Ps") w 0.274(—) r 27 (Ps") = 1.4783643 • 10" 6 r 27 (P S -) 

V 7T' 


(14) 


This formula is based on the result from 


26] for an isolated electron-positron pair (e + ,e ). 


The idea that all annihilation rates with even number of photons are related with each 


other was proposed by Ferrante in 1969 


271] . Later, it was found that this idea is correct 


in application to the Ps - ion and other polyelectrons. Furthermore, the same idea is also 
correct for annihilation rates with odd numbers of the emitted photons. 

The three-photon annihilation rate r 37 (Ps _ ) is ~ 1000 times smaller than the two-photon 
annihilation rate r 27 (Ps _ ). The corresponding analytical expression for r 37 (Ps _ ) is jlfl 


r 37 (Ps ) = n 


4 U 2 - 9 L.5 m 


a b cao i (<5(r+_)) « 2.7185459576 • 10 8 (<f+_) sec 


-l 


(15) 


This formula allows one to evaluate the five-photon annihilation rate in the Ps ion. Indeed, 


by using the formula from 


26] one finds the following result 


/(Y\ 2 * 

r 57 (Ps") « 0.177(—) r 37 (Ps") = 0.955001778 • 10~ 6 r 37 (Ps") (16) 

As follows from Eq. (fl4l) and Eq. (fl6|) the r 4 7 (Ps _ ) and rs 7 (Ps _ ) annihilation rates are sub¬ 
stantially (in cs 10 5 times) less than the corresponding r 27 (Ps _ ) and r 37 (Ps _ ) annihilation 
rates. 

In the three-particle Ps - ion one can also observe the process of one-photon annihilation of 
the electron-positron pair. For an isolated electron-positron pair the one-photon annihilation 
is strictly prohibited. However, a close presence of a third particle drastically changes this 
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situation. The one-photon annihilation rate (or width) Ti 7 is written in the form 29] 


L see 


also 1161 and 


28]) 


647T 2 


r l7 = —— • a 8 • c • flo 1 • (<5 3 2i) = 1065.756921658 • (<5 32 i) sec 

Zj ( 


-1 


where the notation ($321) stands the expectation value of the triple delta-function (in atomic 
units) computed for the ground state of the Ps - ion. This value is the probability of finding 
all three-particles at one spatial point. 

The sum of all partial annihilation rates, i.e. r 27 (Ps _ ), T 37 (Ps _ ), r 47 (Ps _ ), ... is the total 
annihilation rate, or r(Ps _ ). For the Ps - ion we have 


r — r 27 (Ps ) + r 37 (Ps ) + r 47 (Ps ) +... ~ r 27 (Ps ) + r 37 (Ps ) 


= 27ra 4 ca 0 1 


a 


17 19tt a1 


7r 


12 


)J (<y(r+_)) » 100.617460050 • 10 9 (<4_) sec~ l (17) 


As follows from these formulas to determine all partial and total annihilation rates we need to 
know the expectation values of the electron-positron delta-function (<5 + _) and three-particle 
delta-function (<5 32 i) for the ground l 1 *?—state of the Ps - ion. By using the expectation 
values of these delta-functions from Table II we have found the following numerical values of 
all mentioned annihilation rates: r 27 = 2.08048530684 • 10 9 , r 37 = 5.6364151625 • 10 6 , r 47 = 
3.0757152758 ■ 10 3 , T 57 = 5.382786501, T l7 = 3.82491558 • 10" 2 and T = 2.08612172200 • 10 9 
(all values are given in seer 1 ). These values are slightly better than our previous values 


given m 


£• 


Different aspects of the positron annihilation of the t 




in a large number of earlier studies (see, e.g., m 


, Q, 


ie Ps ion were discussed 


m, [29i, 


31] and others). 


IV. PHOTODETACHMENT 


Photodetachment of the Ps - ion is of great interest in applications to astrophysics and 
propagation of radiation in our Galaxy. As is well known the center of our Galaxy contains 
a number of sources of the annihilation 7—quanta with E 1 « 0.511 MeV (see, e.g., (4] and 
references therein). This indicates the presence of objects with very high (local) temperatures 
T > 350 - 400 keV and formation of large numbers of the electron-positron pairs (e _ ,e + ), 
Ps - and Ps+ ions, bi-positronium Ps 2 and other polyelectron species. Photodetachment of 
the Ps~ ion and other polyelectrons leads to very intense absorbtion of the infrared radiation 
in such spatial areas. Photodetachment of the Ps - ion(s) was considered in [5], [32I . Some 
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closely related problems, e.g., elastic electron-positronium scattering and photodetachment 


of the Ps ion by a model Yukawa-type potential, were discussed in 


33] and [34j, respectively. 


There is a well known experimenta paper about photodetachment of the Ps ion 


35]. 


To simplify theoretical analysis and numerical calculations in this study we shall apply 
an effective method of ‘asymptotic photodetachment’, or photodetachment of the Ps - ion at 
very large distances form the ‘geometrical center’ of this three-body system. This method 
was originally proposed by Hans Bethe in 1935 when he considered photodetachment of the 
deuterium nucleus. This method is based on the fact that the wavelength of the ‘acting 
light’ A is much larger than the effective geometrical size of the system a, i.e. we always 
have A>a. Therefore, we can use the wave function of the Ps - ion in its asymptotic form 
which is defined only at very large (or asymptotic) distances between outer-most electron 
e~ and neutral central cluster Ps. Analytical form of the wave function of the Ps - ion can 
be found from the following formula 


T(r) |= Cr-y 1 exp(— 7 ?’) 


(18) 


where C is some numerical constant, Z — Q — N e + 1, where Q is the electric charge of the 
nucleus and N e is the total number of bound electrons. The parameter 7 in this equation 
equals 7 = y/2 A and I\ = x.\ is the (first) ionization potential which corresponds to the 
dissociation of the Ps - ion, i.e. Ps - = Ps + e~. For the Ps - ion one finds in Eq. (lT8lk 
Q — 1, N e = 2, and therefore, Z = 0 and the long range asymptotic of the wave function, 
Eq. flTgjh is represented in the Yukawa-type form | T(r) |= f- • exp(—yr), where C is a 
constant which must provide the best correspondence of Eq. (fT8]l with the highly accurate 
wave function of the Ps - ion at large r. The highly accurate wave function of the Ps~ ion 
is assumed to be known from numerical computations of this ion. 

The photodetachment cross-section a(Ps - ) of the Ps - ion is written in the following form 
(derivation of this formula is discussed in 36]) 


da {H ) = I e / ■ /| 2 do = aa o7^ I e / ' / | 2 do ( 19 ) 

where do = sin ddBdcj) is an elementary volume in spherical coordinates. The normalized 
wave function of the incident state is Cr t~ 1 -exp(—yr), while the wave function of the final 
state is exp(«p e • r). By substituting these exressions into Eq. (lT9l) one finds 

da(Ps-) = aal( ^^ n ^ e )(ef ■ n e ) 2 - I — f ipi(r) sm(p e r)rdr I 2 do 
v 3a; 7 p e Jo 


d 
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= aa n 


, 128irpl 
3u 


r+oo 

(er ■ n e ) 2 - I — / exp(— jr) sm(p e r)dr I 2 do 

Pe do 


( 20 ) 


where p e = p e n e and n e is the unit vector which determines the direction of propagation 
of the emitted photo-electron. Here and below a is the fine structure constant, do is the 
Bohr radius and 7 = \/2I\. The (first) ionization potential I\ of the Ps~ ion is proportional 
to the difference of the total ground state energies of the Ps - ion (see Table 1) and Ps 
two-body system (-0.25 a.u. exactly). The exact expression also contains the factor |, i.e. 
/, = §(0.262005070 .. .-0.25) and 7 = ^§(0.262005070...- 0.25). A few steps of additional 
transformations (see, e.g., [36]) lead to the following final formula (in a.u. u = |p 2 + I 7 2 ) 


dcr(Ps ) = ^ • C 2 aa^( 
9 ' 


Pe 


Pe+T 


r) 3 (k/ x n e 


2 do = 


256vr 2 
9 


• C 2 aa%( 


Pe 


pI + T 


’ sin 2 QdQ ( 21 ) 


where sin 0 = kj x n R and both vectors kf and n P have unit norm. For the total cross-section 


0 


one finds 


/n — \ 5127 T 2 2 / 

-(Ps ) = —— -C «a 0 ( 


Pe 


pi + r 


38.2449007-10 


-18 


r 2 ( Pe ) 

Vr,2 , ^2; 


Pe+T 


( 22 ) 


where C is expressed in atomic units. The formula which allows one to determine the 
constant C takes the from 


N 


C = Da 5 i?exp( 7 j R)T(A,0, R) = D a * Rexp^RjDr 1 

+ exp(—/^.R — jiR) 


2—1 


exp(—ctji? - jiR) (23) 


where D r is the normalization constant for the radial part of the total wave function of the 
ground state of the Ps - ion, while D a is the normalization constant for the angular part of 
this wave function. However, there is an additional normalization constant which is defined 
for the angular part of the total wave function. This constant equals D a = - 7 =. In our 
calculations with the use of 700 basis functions, Eq. (]3|) . we have found that the constant C 
equals C ~ 0.9322567 ■ - 7 =. Now from Eq. ()23]) one finds C = 0.18595831 (compare with [5]) 
and expression for the photodetachment cross-section takes the form 

a(Ps-) = A-(—^-)' = 1.3225275-10- 18 -(—^-) 3 cm 2 (24) 

X Pe + Y K Pe + 7 7 

Analogous numerical computations with the use of 3500 basis functions (see Table III) lead 
to the following value of C ~ 0.1859599866 (at R = 70.5a.it.). The photodetachment cross- 
section is written in the form 

cr(Ps~) = 1.322555 • 10 ~ 18 • ( 0 Pe 0 ) 3 cm 2 (25) 

K Pe +Y 
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The formulas Eqs. (I24j) - (|25|) solve the problem of the photodetachment of the Ps - ion in 
the ‘asymptotic approximation’, which is sufficient for various applications in astrophysics. 
Note that the general approach to the photodetachment of the Ps - ion, i.e. approach which is 
not based on the ‘asymptotic’, long-range approximation of the wave function, has not been 
developed yet. In real applications this ‘asymtotic’ approach can be generalized to include 
more complex cases of the photodetachment of the Ps - ion, including semi-relativistic case 
and photodetachment when in the final state one finds a free electron and secondary photon 
of smaller energy. 


V. CONCLUSION 


The bound state properties of the ground state in the Ps - ion are investigated with the 
use of highly accurate results from recent calculations. At this moment we can conclude 
that the positron annihilation in the Ps- ion is a well studied phenomenon. In the next 
studies we plan to consider some unsolved problems which are currently known for the Ps - 
ion. A special attention will be given to calculation of the lowest-order QED corrections for 
the ground (bound) state of the PsC ion. The lowest-order relativistic and QED corrections 


for the Ps were considered in 


3, 0 


and 


. In particular, in our earlier study 


fl 


we have determined the Bethe logarithm for the ground state in the Ps - ion ( lnK 0 (Ps ~) 
= 3.00502533(5)). The two (singular) expectation values ( r7 3 ) and (rli) for the Ps - ion 
are also known to very good numerical accuracy (see, e.g., 21] and/or Table II). However, 
by using these values and Bethe logarithm we cannot determine the lowest-order QED 
correction for the Ps - ion in the same way as we did for tow-electron atoms and ions (see, 


e.g., 39]). The reason is obvious, since the corresponding ‘small parameter’, i.e. the ratio of 
the electron and positronium masses (m e and M Ps , where Ps is the central cluster), equals |. 
The corresponding series for the lowest-order QED correction converges very slow. Currently, 
it is clear that any actual progress in computation of the lowest-order QED correction for the 
ground state in the Ps - ion can be based on the derivation of the Bethe - Salpeter equation 
for weakly-bound non-Coulomb three-body systems [40] (see also discussion in 23]). 
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TABLE I: The total energies E of the Ps ion (in atomic units). 


K 

E {Ps ; variant A) 

E {Ps ; variant B) 

3500 

-0.26200507 02329801 07770398 027 

-0.26200507 02329801 07770398 256 

3700 

-0.26200507 02329801 07770399 455 

-0.26200507 02329801 07770399 614 

3800 

-0.26200507 02329801 07770400 032 

-0.26200507 02329801 07770400 078 

3840 

-0.26200507 02329801 07770400 250 

-0.26200507 02329801 07770400 279 


TABLE II: The expectation values of some propeties (in atomic units) for the Ps ion. The 
notations ‘+’ and ‘3’ designate positron, while the notations and T’ and ‘2’ stand for electrons. 


<r£> 

(rZl) 

A7 

+ 

S- 

(r=I) 

0.279326542224953365 

0.0360220584545367961 

0.33982102305922030648 

0.15563190565248039742 

(r+~) 

(r —) 

(rl-) 

(r 2 -) 

5.48963325235944993332 

8.54858065509918611146 

48.4189372262379554122 

93.178633847981329005 

(r+~) 

(ri_) 

H-) 

(r 4 -) 

607.295629623278442191 

1265.5804478781441204 

9930.638679796004154 

21054.45338925835809 

([^32^3l] _1 ) 

([r_|—r—] _1 ) 

([^32^31?’2l] _1 ) 

(<5(r + __)) 

0.090935346529989403768 

0.060697690288581955165 

0.0220342380163358569 

3.58891929212-10 -5 

(<S( r +-)> 

(a) 

v+- 

)) 

(a) 

V _ 

0.020733198005219 

-0.500000000057701 

1.7099675635500T0 -4 

0.4999999981459 

TeN 

7*ee 

(/) 

CO 

•-s 

CO 

to 

coco 

0.5919817011489022333 

0.01976963281713200176 

0.0509332587787341171 

-0.1234320911052318 

Hv 2 _> 

(~W) 

(Vr • V 2 ) 

(v+ • v_) 

0.0666192945358900085 

0.1287664811612000907 

-4.47210791057992633T0 -3 

-0.1287664811612000907 

( r +- ' r__ ) 

CO 1 

1 + 

CO 1 

1 1 

CO 1 

1 + 

46.5893169239906645026 

-0.25348417470871 

0.01131050073132 

7.056875445764-10 -3 


( a )The expected cusp values (in a.u .) for the Ps ion are z/_|_= —0.5 and v _= 0.5 (exactly). 
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TABLE III: Numerical constant C from Eqs.(22) (23) determined for different electron-positronium 
distances R (in atomic units) for the ground state of the Ps _ ion. The constant A is defined in 
Eq.(24). 


R 

C 

c 

s/Sn 

A 

65.5 

0.93224641707 

0.1859562558 

1.322502E-18 

66.0 

0.93224854366 

0.1859566800 

1.322508E-18 

66.5 

0.93225060763 

0.1859570917 

1.322514E-18 

67.0 

0.93225261127 

0.1859574913 

1.322519E-18 

67.5 

0.93225455679 

0.1859578794 

1.322525E-18 

68.0 

0.93225644627 

0.1859582563 

1.322530E-18 

68.5 

0.93225828175 

0.1859586224 

1.322535E-18 

69.0 

0.93226006515 

0.1859589782 

1.322540E-18 

69.5 

0.93226179830 

0.1859593239 

1.322545E-18 

70.0 

0.93226348299 

0.1859596599 

1.322550E-18 

70.5 

0.93226512090 

0.1859599866 

1.322555E-18 

71.0 

0.93226671366 

0.1859603043 

1.322559E-18 
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